ANAi VI 1C UBVtiOPttt NT Of THE FLIGHT PATH 

equations lor maximum range by the 

METHODS OF VARIATIONAL CALCULUS 

NHtMUM «M» >f4r 



1 II: ID LNANT COMMANDER CHARLES W. FRITZ, USN 
LIEUTENANT JOHN M. TIERNEY, U8N 




ANAL me DEVELOPMENT OF THE FLIGHT PATH 
EQUATIONS FOR MAXIMUM RANGE BY THE METHODS 
OF VARIATIONAL CALCULUS 



by 

Lieutenant Commander Charles W* Fritz, USN 
Lieutenant John M 0 Tierney, USN 

Aeronautical Engineering Report Number 305 
Kay 20, 1955 



Submitted in partial fulfillment of the requirements for the 
degree of Master of Science in Engineering from Princeton University, 1955o 



ACKNOWLEDGEMENT 



The authors wish to acknowledge the encouragement, 
guidance and assistance extended to them by Professor D,0* 
Dommasch of the Department of Aeronautical Engineering of 
Princeton University, They are indebted to him for the 
concept of this thesis which is based upon some of his original 
work in the adaptation of variational calculus to aircraft 
performance problems© 







TABLE OF CONTENTS 



PART I 
PART II 
PART III 
PART IV 

PART V 



Introduction 
Development of Equations 
The Bounded Case 

Evaluation of the Derivatives and 
Formulation of the Final Equations 

Method of Solution of the Equations for 
an Extremal Fligth Schedule 



CONTUSIONS 

APPENDIX A Reduction of Euler Equations 
to Classical Range Parameters 

APPENDIX B Consideration of Mach 
Number Effects 



BIELIOGRAPHY 



page 1 
page 6 
page 25 

page 31 

page 38 
page h5 

page h6 

page h8 
page 52 



AELREV1 AXIOMS A.'H S^TQjlS TO I 2 USED IN TUf RETORT 



Wing aspect ratio 

Aircraft total dra D coefficient 

Induced drag coefficient 
Drag coefficient at zero lift 
Thrust specific fuel conrurrrotion 

, Centrifugal force 
Total aircraft drag 
Total differential operator 
Partial differential operator 
Induced drag 

Drag at zero lift 

Airplane efficiency factor 
Total effective thrust 

Thrust due to momentum change of air 
Effective thrust available 

Gravitational acceleration constant (32.17 ft/sec^) 
Altitude 

Aircraft lift > 

Mass of air and fuel 
Mach number 

Dynamic pressure (q - 1 p V 2 ) 

T7 

Wing surface area 
Time 

Engine exhaust velocity relative to the origin 
Exhaust velocity relative to the body (V - u) 



Velocity 



Speed of sound 

Velocity derivative with respect to time 
Total weight of Lod} or aircraft 
'..'eight of fuel expanded from point 1 to 2 
Horizontal distance 

Flight path angle relative to the x-axis 
Flight path angle derivative with respect to time 
Air density 

Angular velocity in pitch 



Aircraft total lift coefficient 
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and altitude, such as pt. A on fig. 1, and flovm a given distance Xj_. 

Tan schedule of altitude and soeed would be determined l ; the inter- 
relations of the variables of the problem bounded by the capabilities of 
the particular aircraft. The problem could just as well 1 e started 
fom the instant of take-off, sa; pt. E, and thus cover the entire flight 
fron the take-off until over xq. 

Commencing a maximum r^ng 5 or minimum fuel schedule from 
take-off is a completely new concept. Classical range solutions do 
not attempt to solve for an optimum flight path from point of take-off 

to optimum cruising altitude, but ra trier specify, for turbo-jets,*^ 

* / 

a .nebulous "minimum fuel climl " to the highest altitude rt which 

speed for (Cq /c [)) ffiX is equal to V nnx . The aircraft then cruiser 

at tnet altitude, slowly climling as the weight is decreased, thus 

Vl . 

increasin’, tnc altitude at which the speed for maximum Cj /Cm is 
equal to the maximum speed. 

This "minimum fuel climb" will, if properl; determined, 
achieve a certain altitude with a minimum fuel expenditure, but this 
is not the primary concern. The problem is to travel a distance xi 

9 

with minimum fuel consumed. If the variables of the problem call for 
a climb to a specified altitude, then that climb should be executed to 
give maximum distance for a given fuel expenditure to get te that 
altitude. As long as the flight path stays within the limitations of 
the aircraft, however, there should be a single path with no dis- 
continuities for minimum fuel path for a given distance. Thus 
looking at fig. 2, it can be seen fiat to travel from pt. a to some 
pt. £ a distance Xq avray there must be some minimum path ac which usqs 



less fuel than the path abc wtiich is composed of two parts, even 
though aL and be are sopai'ate paths. 

Classical range equations, moreover, are not generally 
applicalle or accurate for present day high performance aircraft. 
Classical range performance is based on assuming that the flight 
path angle ( 4 ) is small such that its cosine ic considered one, 

(this assumes that L ■ W) and that flight path accelerations along 
the path nere ignored (i.e. Thrust a hrag). From a review of the 
literature on the subject, it was determined that everything written 
to date on the range of aircraft nos been based on t :ese sjr-r. lifting 
but very restrictive assumpti ons. With modern high performance air- 
craft, flight path angles up to and including 90 degrees are possible 
and flight path accelerations can easily exceed 1.0 g. In fact, in 
considering the gener''l problem of powered flight (rockets, missels, 

VTO, etc.), the convenient and conventional condition of steady iev^l 
flight is almost the exception instead of the. rule. 

It is the intent of this work to indicate methods based on 
the principals of variational calculus for tne determination of schedules 
for minimum fuel or maximum range flights for all types of powered 
flights. It is intended, more specifically , to develop equations 
for the minimum fuel paths of turbo-jet type aircraft. This differs 
from th"' classic problem in that the schedule begins at take-off and 
minirizes the fuel expenditure ffo.r that point to a distance Xy, taking 
into consideration possible large values of flight path angle and 
considerable flight path accelerations. In other words, the problem 
considered nere is the broad ;;enerr.l problem. It uses accurate, 
realistic conditions an 3 minimizes fuel consumption for the entire 



flight path from 11 wheels up ,: to "over target". 

It is hardly necessary to discuss the imoortance of reduced 
fuel consumption since its advantages are so readily apparent, but 
to help keep these advantages in mind, a few of then will be listed: 

(1) Increased useful payload (bombs, cargo, passengers, etc) 

(2) Greater range or combat radius with the same' payload 

( 3 ) More economical operation overall 

(10 Complete specification of a flight schedule which could 
be incorporated into an automatic control for a manned 
or pilotless aircraft. 

The adaptability of these and other advantages to all phases of 
aviation is al30 readily apparent. The long range bomber, the short 
range fighter, the commercial cargo and passenger carrier, the pilotless 
aircraft, and practically all other vehicles traveling through the 
air could adopt flight schedules determined in this manner to increase 
their effectiveness and their efficiency. 

A specific application for which this work could be im- 
mediately useful is the determination of entry and escape paths for a 
supersonic bombing mission.. Those parts of the entry and escape 
path not determined by other tactical considerations would most prob- 
ably be flown at full power including afterburner and following a 
path to give best fuel economy at maximum power. The bounded thrust 
path equations developed lata* in this thesis apply directly to this 
prollem. 

In this and many other applications, Kach effects must be 
considered. For simplicity, however, Mach effects have been omitted 
from the derivations that follow, but have been treated separately in 
Apoendix B* 



PART TWO: U2V.iL OP! ENT OP % A HONS 



The range of an aircraft is, as we have hypothesised, a 
function of its flight path. Its flight path, in turn, depends on 
varial les such as amount of fuel, velocity along path, flight path 
anjle, accelerations along and normal to the path, and thrust spec- 
ific fuel consumption. These varial les are, in turn, functions of 
each other, drag, thrust, lift, weight, and air density. Relating 
these variables among sufficient equations to solve for a minimum 
fuel path is the problem of this thesis 0 

A) EQUATIONS OF MOTION 

The equations of motion pertinent to the problem are the 

force balance equations of a body moving in a vertical plane. This 

motion will include accelerations' along and normal to the flight path, 

as in maneuvering flight, since these accelerations will, a priori, 

affect the loads and thrust required, and therefore affect the range 

(or fuel consumption). The forces acting on a body moving along an 

arbitrary path in two dimensional space are shown in fig. 3. 

If this body moving space is a conventional aircraft, 

% 

the forces F]_, F2, and F3 can be expressed as: 

Fq The total effective thrust. (F) 

F2 The lift of the moving body. (L) 

F3 The total drag force. (Body drag, D, and inertia, n dV ) 

. ^ 

Ey equating forces parallel and perpendicular to the flight path to 
zero, the following equations are obtained: 

(1) Forces parallel to flight path 
£ F (( » 0 ■ Fq - F3 - WsinV 



( 1 ) 
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DuTANce; ; X 

F.g 3 

poRCfj ACTING. on Body 






or for an aircraft: 

F*-D - d (mV ) - Wsin*' ♦ u *2_ - 0 

3t^“ dt 

where F, D, W, , u, and ra were defined previously and: 

d(mV) ■ n dV ♦ V dr. 
dt dt dt 

But 

dm ■ -Oj>F 
dt ~g~ 

since fuel is the only variable mass involved* 

Also dV = V 
dt 

Now 

£F ■ F*-D - mV + VGjiF - Wsin^ - uC^F • 0 

1 “ ~ 



( 2 ) 



( 3 ) 



(U) 



Considering the definition of thrust, it can be seen that the total 



effective thrust acting on the aircraft would be equal to that im- 
parted by the air through the engine as it experiences a change of 
momentum plus that imparted by the fuel mass leaving in the exhaust 
due to its momentum change* The thrust due to the leaving fuel is a 
function of the difference between aircraft velocity and exhaust 
velocity and is expressed by: 



Thus 



(u - V) Of F 
g 

F - F* - (u - V) Ct F 

“IT 

= F* + u e Cy F 

"IT 



( 5 ) 



The equation now becomes 



F - D - W V - Wsin i - 0 (6) 

g 

(2) Forces perpendicular fco the flight path 

- 0 ■ F 2 - CF - W cos-J (7) 

? 



or for 


an aircraft 






L - CF - Vcos^ * 0 


(8) 


But 


C F * ra V u) 




and 








u> » d{ s 






dt 


» 


then 


# 






CF - w vi 


(9) 



g 

And therefore 

t 

EF^ - L - W V? - Wcos/ - 0 (10) 

I 

Equations (6) and (10) above are the general equations of 
motion pertinent to flight in the vertical plane. These will be modi- 
fied later for use in the range problem to the extent that derivatives 
with respect to time will be replaced by parametric derivatives with 
respect to either weight or distance. 

B) EQUATIONS OF MOTION FOR RANGE 

Eefore attempting to maximize the general equations of motion 
we will further modify them to be more applicable to the problem of, 

(1) Maximizing the range, x, for a given fuel expenditure,^ 1 ^* or 

(2) Minimizing the fuel expenditure, jWgfor a given range, x. 
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For maximizing range, x, it is convenient to '•'} J ani>p variables 
as follows : 



define 


d( )/d! ( )• 








Then 


V - dV - N ci; - 


y __ 








dt <i; dt 


dt 






But 


J - -Gp F 
dt 








Therefore 


V - -7» Sp F 






(11) 


^lSO 


^ * dtf = ch( 


= 1 


a-; 






dt d ' 


dt 


dt 






^ - - t<» C T F 






(12) 


Substituting 


(11) and (12) into (6) in 1 


(10) gives: 




L + 


1 V if* C T ^ - :;cos)( - 0 




(13a) 



F - T) 4 U V’ C? F - ,T sin^ *» 0 
K 

The integral equation relating x and in simplest form i3: 



(13b) 





X - 


/dx d f 






; ** 


Now 


rbr b 


dt 




U. . 


dt d- 


But 


: = 


Vcos'/ 




dt 




find 


fron before 


i 




dt - 


- 1 




d. r 


Oy F 


so 


dx 0 


- Vcos ^ 




!z: 


Gt f 


And 


therefore 






x » 


/ - Vcos 






j C T F 


From 


J 

i (13b) 






TU 

r “ 


0 >. Ms in 



( 1 ?) 



1 + _ '"r 

R 



FcosX - I 

•^fvt vr 
J/ |{ 



( 16 ) 



(17) 



Substituting into (15) 



x 




dW 



or 



x 



Vcos/ (1+ Wl'Op) dW 
Cp (EH-’Wsino ) 



( 18 ) 



(2) For minimizing the fuel used, jVr'2, the changes of 



variables 


found useful are: 


Defining 


IQ; ( )' 

dx ’ 


Then 


V - dV - dV dx + V'Vcostf 
eft dx dt 


And 


$ * d# * dj( dx' » ^’Vcos^ 




dt dx dt 



(19) 

( 20 ) 



Substituting into (6) and (10) gives: 

cos - Wco3j(= 0 
g 

F - 15- V/ V' Vcostf - Wsin^ - 0 

& 

The integral equation relating x and W in simplest form for minimizing 




As before 




and 



dtf 

dx 



or f 

Vcosl^ 




Op F 
Vcoslf 



dx 



(22) 



W 



Solving (21b) for F and substituting gives 

W - j - Or (D + ^V' Vcost( ♦ WsinjQdx 



/ 



V costf 



or 



W - 



V- [- OS r D * WT< ♦ '.'tan/ ] , 

Now, although equations (18) and (23) relate the variables 






+ wvi + 

g 



(23) 



of the oroblem, the; do not specify the flight path. They are equations 
for the general motion of a body in the vertical plane, expressed in 
such a way as to give the relationship between range and fuel ex- 
penditure for any flight path. The variables under the integrals 
are dependent on the path, yet the equations as such do not specify 
a path. To determine equations for the particular path which will 
give the greatest range or least fuel expenditure ^t is necessary 
to look to the mathematics of maxima and minima. Presenting the 
problem in more mathematical terms, the unique path desired is that 
path giving the extremal value of an integral. The problem of finding 
an estremal value belongs to that branch of mathematics known as "The 
Calculus of Variations”. A brief discussion of this subject follows. 

C) THE EULER EQUATIONS OF THE CALCULUS OF VARIATIONS 

Let I be a quantity defined as an integral function of 
several variables as follows: 



where ( )' indicates the first derivative with- respect to x. The 
Calculus of Variations is concerned with finding the extremal values 
of the quantity so defined,,, The necessary and sufficient conditions 




LL± 




F '% 4 

FAMILY OF PATHS FROM 
POINT "0" JOT 
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V (x) “ 1' (x) + *■ n 1 ?' ( x ) 

With the conditions: 

\i*o) - f(Xo)) 

_ f i.e. V7(x) 

V x i> a Y < x i)J ' 



0 at the end points 



(26) 



Where € n is a small constant quantity and r^(x) is an arlitrary 
function of x that is continuous and has continuous first and second 
derivatives in the range of integration. 



3^(x) therefore defines all possible paths of integration 

between X. 0 ,Y 0 and X^Y^ so that (2li) may be written: 

_ _ 

I n “ / 0 fx,Y n (x),Y n » (x)] dx 

or replacing Y n and Y n ' from (25) and (26) 

'X„ Y 



n “ f 0 fj[,Y(x) +€ n ^(x),Y* (x) + € n r^'(x)] dx 
J *cX 



( 2 ?) 



It is - clearly obvious that I ^differs from Ig only through the value 
of € and it may be concluded that: 

1 n “ 1 ( € ) * * (O only 

and since € is not a function of x, the ordinary maxima ahd minima 

l 

rule may be applied, i 0 e c : 



dff 

d€ 



Hi » 0 



*-*•0 



for an extremum 



and, furthermore, the differentiation may be made under the integral, thus: 



dG ' 

dT d* 



( 28 ) 



dG - ^_G , d f X(x) + c r r>(x)] 

de <J|Y(x) ♦ £»»7(x)J <*<=■ 



From (2?) 



Since 



1 ' 



♦ ^ G 

dlY'(x) + e n 7'(*>J 



clfY 1 (x) + u yy (x)1 ♦ dx 

fit d£ 



dx - dY * dT? ■ dl' * dv' * 0 
de d£ de dt d£ 



then 



or 



™ “ d^i . 7j(x) + dG . t?' 

d€ d(_Y(x) + € n >j(x)J ' d\v (x) + € n (x)J ‘ 

/ f c)G n + dG "W* 1 

J L^ y+ M ) +6 n V) J 



6o 



(x) 



dl 



dx * O 
(29) 



and 



/■* i»X /'^•J 

dl » / £G 77 dx + / oG Y7* 

^ y d(y> £ n ^) ' y ^(y + € n Y|') ( 

A Jo **> Y . 



dx = 0 



From whence 



*,X 



dG V7 dx • - i dG y'dbc 
c3Y n 1 / At . 



^Y_ 



%,y, M, 



Integrating the right side by parts and remembering *that ^ (X D , Yq) 

7 M ' 9 , 



c)G w-'dx ** - 



dG_7 

^n‘ ( 



w 

"cEx 



(30) 



7 



dx 



*.J. ✓x.x 

Substituting (30) into (29) gives: 
-*>X 



dl 

dl 



~ _d ^G >i dx » o 



dx 






( 31 ) 



*o,Y. 



Evaluating (31) for € * 0 to find the extremum: 

x,0C 






**X 
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But this is true onlj if 

|>G - d (AG \|v- 0 
[_3x d^VOTj 1 

throughout the range of integration. Since y is completely arlitrary 
at every point except the end points, it necessarily follows that for 
an extremum: 

dG - _d fdG ) m 0 (32) 

c)y dxl^y 



This equation is called the Euler equation. It must be satisfied if 
the integral is to have an extremum value and it is through this 
satisfaction that the path, Y(X), is found. 

Proceeding to the more general problem of several varialles, 

the integral must satisgj the definitions as before and: 

x Y Z 
Vi o 



G [X,Y,Z, , Y',Z', ] dx 



r . y 7 ' 

Proceeding as before ^ we arrive by an anaiagous developement at the 
Euler equations: 



AG - d f AG \ * 0 
3y 

<)G - d {AG \ 

dz ^ l3z«/ = 0 

<)G - d &G - 0 etc. 

3T7 *E“ZT) 



A 

► 

J 



(33) 



These equations (33) must be satisfied simultaneously for an extremal 
value of the integral, such satisfaction therefore defining the desired 
path as in the two dimensional case. 

It should be remembered here that an extremal value for the 
integral may be either a maximum or a minimum value. Thus in deter- 



V 



mining a solution for equations (33) > we nust determine ly other . t 
means whether a particular solution will yield a maximum path or a 
minimum path* For most types of engineering work this will be in- 
tuitively apparent* 

D) T4E EULER EQUATIONS FOR KAXIKUK RANGE 

a.) Minimum Fuel Equations (i.e«, x independent variable) 

Equation (22) is of the form 

- X W 2 - 1(c) “ / Gdx 

•A 

where 1(c) = -W the quantity (of fuel) to be minimized and G = Op F 

Vcostf 

the integrand function* 

Now: 

G » Q(Ct,V,*,F) 

F - F(D,V',W,V,Z() 

D = D(L,V) 

L =* L(W,V,rf,rf») 

thus 

G - G(W,V,*,Qr,tf»,V') 

Now, C.p, the tnrust specific fuel consumption, is a variable, and a ve.ry 

important variable where range is concerned, but treating it as such 
in the derivation of Euler equations gives rise to extremely com- 
plicated relations and to derivatives on Op that are very difficult 

to evaluate. We have therefore elected to treat Gp as a stepwise 

variable constant* 

G - G(W,V,Tf,V',tf' ) 



Thus 












I 1 

1 



And the Euler equations of this function become: 



1.) G - d fdd \ » 0 

“3v dx vTvv 

2* ) c) G - d fd(} \ - 0 

■a? 

3 .) A! - / 3 g >| = o 

dW dx(Jw r / 



1 .) Derivation of the first Euler 
G — d f dG \ • 0 

TT $U v 7 / 

Now from ( 22 ) 



G - C T F 
Vcos^ 



Then 



And 



<^G 

“ 5 v 



Vcos^ Op - Op F cos^ 
— -y JV 

(VcOST^ 



g r^v - 1 \ 
L F vj 



<) G » °T /^F 

<jV Vcos^^v y 



_d ^G 
dx dV 



W <* df 



c)F Y y« cos^ - Vri'sin ?6 



(Vcos^)^ 






.t 



Using an abbreviated notation: ° j = Y„ 

o 2 L 



_d ^ dG V C p 

ax Uvy p 



Fy» - F yt ( V* cosl( — 
F“ FYcosrZ 



•sin^)| 



or 



Gyi m G |~Fyi 1 — V'Fyi + Fy» ^’tanK^ 



V F 



V * I 

F J 
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G Fyi 
F 



Fyi ' - V' + /' tan/ 

PyT T 



Now! 



Gy - Gyi 1 - 0 * G 



fFv 


- G Fyi 


L F xJ 


— 



~F yt » - V* + /' tan/] 

L*^- * J 



Multiplying through by F gives us the first Euler 

f G 

F y - F — Fyt * + V* Fyt — Fyi /'tan/ ■ 0 

v T* ' 

2«) Derivation of the Second Euler equation 
d G • d d G m f) 

17 53? 

as before* 



(3t») 



and 



then 



G » Op F 
Vcos/ 



c)G • Vcos/ dMl + Op F Vsin/ 
^ ( Vcos h j 2 



F^ ♦ tan/J 



c>G - Cr ft ? A 
d/* Vcos^^ 7 / 



dl&G 

c£x\<5/ 


r) “ Gp 


— ~i 

~Vcos/ gr' - F^» (V cos/ - v/'sin/) 






_ (Vcos/)^ 



d dG 

dx np ■ 



- G 


1 

c* 

«• 

l 

r> 


r« - /» tan^ 




F F \ ' 


f / 



G JV 

F~ ; 



F^.« - V + 
-1- — 



/'tan/ 



lV 



Nov; forming the Eulers 



o, - cy - o 



tan/ 



GJV I v 

F [ F s‘ 




Multiplying through by F/G leaves the second equation* 

+ Ftan/ - F^,' ♦ V* F^« - F*. tf'tantf - 0 (35) 

3.) Derivation of Third Euler 

dG - _d Ag ^ ° 0 

dW cbcl^W 7 / 

G - Gp F 

VcosV 

<)G * Op /AfI 

"oW / 

dG “ 1W» - 1 

■Jv* *3w T 

_d_/^G_\“ 0 
dx law 1 / 




(36) 



Therefore F^- ■ 0 and there is no Third Euler in this case. 

b) Maximum Range Equations (i.e. W is independent variable) 

Note* d( ) ■ ( )' ’ 

cW 



From equation (15) 

G * - Vicos'/ 

Op F 



1. ) Derivation of the First Euler Equation 



dG 

IT 
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Forrcin, 



And 



<^Q - - Vcos/ | - <)f/<)V' 

Tv 1 " or [ (f )2 



Vco 3T^ Fy, 



op F" 

d (6 G \ “ C-n F 2 (V* cosi Fy, 



Vcos^Fy, • - V /'sin ~i Fy«J 



Vcosj'Fy, 2F F* C T ^ 

\ «>r F2)2 y 

V* Fy, cos}{ + Voosi Fy, * - F v , V D'sin/- 



Gy.' 



Gp f 2 

» 2Vcos/ Fy, F* 

Or f^ 

■ G f V* Fy, 

[ v-p~ 



Oj, F 2 Op F 2 



V* Fy, + Fy, ' 



- ^'tan^ - 2F V , F» 



J 



G rFv«lf2F» + tf'tan^ - Fy, » - V» 




V» 



■] 



l the Euler 
Gy - Gyi 



G 



[JL - Fy " 


- G Fy, 


2F« + Y»lW - Fy«» - V* ' 


[v T ’J 


~ 


jr f v » -Tj 



* 0 




2F V , F« 
F 



F y, 



^'tantf + F 



V * 1 



t 

+ V* Fy, - 0 

y 



( 37 ) 



2.) Derivation of the Second Kuler Squation 



dG 

IT 

& G 

IT 



d fjG \ - 0 
aw vST 7 / 

V | » F sin"/ - cos / pF/jJ ) | 
tan"/ + 



Op 



8C» - Vcos/ r. » 

TF Ct 



(F)2 



-VF^< cos / 



— /^JL^ 
dW UW 



- JL_ [fy cos/ V* + Vcos/Fj, 4 - V F x , K'sin/ 

°r f2 f5 f2 



<jG 
dW 






- V F^i cos i (2F F') 



a,/ - -G [XL 5* + 

l v F F 



F^tan/ - 2F^i F' 
F ~F2 



Now 



G V - V 



tan A F./ +G|V F <( + F/ - ^F yJ tanjf 
F ’ [fF “F F 



• 2F 2r 'F' 

T" 



0 



From wliich: 



F tan/ + Fy - \M F^ - F^' + /'F ( tan/ + 2F ^F* - 0 (38) 

V V T~ 

Thus equations (3h), (35)# (36), (37), and (38) represent extremal 
paths in terms of flight path angle (^), velocity (V), horizontal 
distance (x) and pounds of fuel (W) • The first three equations re- 
present extremal paths for minimum (or maximum) fuel consumption for 



traveling a given distance x, and the last two equations represent an 
extremal path for maximum (or minimum) distance x for a given fuel 



expenditure W* 



PART III 



2 



THE FOUNDED CASS 

The equations derived in section D of Part II apply to a 
more or less general propelled body with unlimited thrust and thus 
unlimited velocity and flight path angle. They do not consider any 
boundaries on the thrust available to follow the extremal path that 
they designate. This extremal path may very call for altitudes, 

velocities, and accelerations above the physical ability of a parti- 
cular airframe - engine combination for all or part of the path. 

' ,< 

Those equations are, therefore, the absolute extremums for an ideal 
body without limitations imposed on any of its parameters. They 
define a path that is the best possible and should be followed when- 
ever that path is within the operating capabilites of the particular 
propelled body. 

The primary limitation imposed by physical factors in any 
particular case is that of thrust. Thrust varies from zero to a 
finite maximum value that is a function of the individual thrust 
producer and other variables. For a given .type power plant the 
maximum effective thrust is usually a function of the variables of 
the problemo For examples The thrust of a rocket is a constantj the 
thrust of a ram jet varies with velocity and altitude] the thrust of a 
turbojet varies with velocity and altitude 0 Therefore, for any given 
airframe - engine combination there is a limiting envelope or boundary" 
of thrust which limits the path variables such as velocity-, flight 
path angle, altitude, etc. Whenever the general path equations of 



*L~\ 
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Part II demand a thrust that exceeds this boundary, the path becomes 
unrealistic and of academic interest only. It then becomes necessary 
to compute a new path from this point that will be an extremal based 
upon a maximum continuous value of thrust available. Euler equations 
will now be derived to define an extremal path 'for the case of bpunded 
thrust. 



A.) The Bounded Case for Minimum Fuel for Turbojets 

As was said before, the thrust of a turbojet varies with 
altitude and velocity, i.e. : 

F * f (V,h) where f pertains to bounded thrust 

Eut since altitude is a function of path angle, we may also 

t 

say: 

F - f (v,i) 

From before: 

O - 2J., 

V cost 



But now 

g - or f 

Vcos'J (39) 

This function is independent of the characteristics of the airframe - 
engine combination and must be modified to include them 0 From (21b): 

V cos~{ m g(f - D - Wsin?() 



Substituting this into (39) above 

G - Or f W V 1 

g(f.D-Wsin^) 



(to) 



r 
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1«) Derivation of the First Euler Equation 



c>G — d / <) G ^ * 0 

“ 5 T »“ 57 *y • 



<)G 

av 



F bt/bv - Af/AV - AD/^V 1 
L f f - D-WsinTs J 



(ia) 



But 



Then 



Now 



G_ - Qr £ w 
-V' g(f - D - WsinH( ) 



W ' 



- G 

V * 



AG - W » 

“ 3 v r — 



d 

dx 







- vw - W 1 V * 1 
( V «) 2 



(la) 



Forming the Euler (la) - (ii2 ) ° 0 



AG - d ( Ag ^ - W' [* fy - f V - Oy 
av dx V <$'"/ |f (f-D-WsinVj 



ytyt' _ W »y»i , o 

( V *) 2 



From which 



vf” - V' f *7 - -Pv ~ 1 - v” r o 
w * [f~ (f-D-Wsin K)J ”j 



(143) 



iG . 



2*) Derivation of Second Euler Equation 

W-r^. Af/£^ - <)D/)if - Wcos^ 

[. f ( f - D - WsinV ) 



( I 0 ») 



= w»r i 

bV [(f-D-Wsini) I 

• "* « » 

J ( W' f(f-Wfein«V - Dji (f’-D'-W’sin# - Wy^cosjOl 

dx\bi'J L (f-O-Wsin^)il (f-D-Wsintf) 2 J 



+ W" f *V 1 

[ rr^^iTHTi J 



(145) 



Forming the Euler (Uh) - (Ii5) * 0 



W" 



V j 

(i'-S-VsinO 



•■[? 



W f * - - Wcot-rf) 

rr - r^TTsinT) 



+ w 



, 1 / , / , . - Dv/ (f'-G’-V'sin* - W*costf) 

U’-h- «sir^ J 



= 0 



(f-’V.ssirrtO^ 



From v;hlch 

W“ D^' - _fv (f-D-Wsintf) + f^ - 0* - Wcostf + fy* 
f 

- (f '-D’-W’sin^ -W^'cos^ ) . q 

(F-TU&'i n*J (b6) 

Partial derivatives of drag with respect to V* and thrust with 
respect to t£' are zero since drag and thrust are not functions of 
V' an* i>' resoectively. 



E) The Pounded Case for haximum Range (i.e. with W as independent 
variable ) 

As before 

G = Vcosl{ 

Op F 

and 

F - f(V,h) =■ f(V,^) only 

50 

G - Vcos { ' (Itf) 



Introducing the airframe characteristeics s 



From (ija) 



- Wcos/ ■ 0 

V { ' Op f 

g 

(1*7) leaves* 

&A . (1*8) 

g 

t 

1.) Derivation of First Euler Equation 

()G - L + VLy . V lvf + 2 VY' (1*9) 

”"3v" C T Wf C T Wf GpWf2 g 

c)G ■ VLy. (thrust is not a function of V') 

Tv 

d /■ do \ - ”*•* * V'Lyi _ V Lv> (wr * t) leM 

~V~w) ~"6i v:"r ‘ " 0 r(wf)i! 

Forming the Euler (1*9) - (50) * 0 



L ♦ W V /*C T F 

i 

cos'/ • L + 
W 



Substituting this into 

g * vl — + 

OpW f 






L + VLy - VL fy + 2 fV “ VLy i 1 

~ & 



V'Ly, 



+ VLy. Wf' + f 
Wf L 



]■ 0 

(51) 



2.) Second Euler Equation 



^G “ 


V 


'l y - Lfv 1 


Of 


(52) 


IT 


OpWf 


L f J 




j£- “ 
IT. 


C^Wf 


+ _v£ 

g 







✓ 
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d [A C>\ - VV ♦ V* L t * . VLj,[wf’ ♦ f] ♦ 2W 

*~(TrJ Opv/f c^wf or(wf)^ e 

(53) 

Forming the Euler (5 2 ) - (53) " 0 



L, - 



L f* - L#*' - V' L,' + Ly*[wf' + f ] - 20jW f V' 
T" T Wf “ 



(5U) 






PART IV s EVALUATION OF THE DERIVATIVES 
AND FORMULATION OF THE FINAL EQUATIONS 

The Euler equations as developed in Parts II and III contain 
derivatives on several of the variables. To determine the desired 
path, these derivatives must be evaluated and included or discarded 
depending upon their effect. In this thesis only the derivatives 
for the minimum fuel path, both bounded and unbounded, will be deter- 
mined. The evaluation of the derivatives for the maximum range path 
(W independant variable) is analogous but the expressions are much more 
unwieldy. 

-if 

Considering first the derivatives to be evaluated for the first 
Euler in the unbounded case, we have: (l)Fy, (2)Fyt, (3)Fyi'. From 
equation (21b) 

F - D + W V» Vcos i + Wsintf (55) 

I 

Taking the desired derivatives.: 



F y - D v + W V* costf 






(56) 


Fy|® Dy, ^ _W Vcos’tf 




W V'cos^ - i' W Vsin^ 


(57) 


Ftm ' - D„, ' + W'Vcos i 


+ 


% 


g 




g g 


* 


■ D nl ' + W Vcos/j 


1* W» 


+ V' - Tf'tantfl 


(58) 


T 1 


L w 


” J 





Consider now the derivatives to be evaluated for the Euler in K 



for the unbounded case;, the following appear: (l)F^, (2)F^» , 
(3)F*i' # From (55) 

F - D + W V V'costf + Wsin t 
g 



Whence : 



? i " D * - JL v V'sintf ♦ Wcos ^ 




(59) 


g 






" D *' 




(60) 


F^ • • ■ D ^ ' 




(61) 



These, then, are the desired derivatives as functions of drag or of 
drag derivatives. The drag derivatives will be determined after a 
look at the derivatives appearing in the bounded equations* From 
equations (h3) and (1|6) for the bounded case we have the following 
derivatives* 

(1) W« 

(2) W» 

(3) f v 

<W • ' 

(5) f 

Now W' * Op f W V’ 

g(f-D-WsinX) ' (&) 

From which 

w „ Op f W V' f ft + W i + V" - (f - D' - W*sin^ - Wtf'cosJO] 
g(f-D-Wsin#)[ £ W V* (f-D-Wsinfl) J 

(63) 

Since the thrust, f, used in the bounded equations is an actual 
thrust, (i.e* thrust available), engine data for the specific 
engine - airframe combination will, supply the information necessary 
for evaluating the derivatives on thrust. 

i 

f ■ f(V,h) A known relationship for a given aircraft 

Then 

f v 31 bf/bv' This can be taken directly from engine curves (6i*a) 



33 



Now 

f, - df c>h 

‘ ■51TTT 



(61(b) 



and df/<ih is available from engine curves while b is a simple 
geometric relationship© 

f« - <^f V' + hi V' (6Uc) 

Tf TT 



To evaluate the derivatives on drag it is most convenient to 
express drag in the following form: 

D - D 0 + D i 

* 1 /)S V 2 C D + 1 flS'# .a 2 

tv ”/ vsrr 



Or, substituting for its equivalent in terms of lift from equation 
(21a): 

D - 1 />S V 2 C Dq + 2 (Wcos)Q 2 / 1 ♦ V 2 f \ 2 ^ (65) 

° TT eAR yO S V 2 V g / 

In taking derivatives on the above equation for drag, the subsonic 
region only will be considered. No Mach effects will be present 

then and both Cn and e will be treated as. constants. An extension 

u o 

of this analysis to include Mach effects will be found in Appendix B* 
Thus 

DV - />SVCp n . u (Wcos*) 2 fl -/ Y 'V 2 . ) 2 ~ 

1 HeAR V3 |_ ^ g / 

= 2 Do - 2 Di (l - j (66) 





D v » -0 (67) 

Dv‘"--0 (68) 



A 



(69) 



* -It V) cos 2) sin^ f 1 ♦ y2 5* j2 



"Ti elRjoS 
-2 D t tan'tf 



V ^ 



ty 



h (WcosT()2 (x + ;;g — } f V2 \ 

.^eAR/oS V 8 / 



,-ff eARyoS 
2 Dj 

V*' 



(70) 



V 



Wcvti ( x «. V£j 

iTeJR pS g 



2W» - 2^'tan^ ♦ (V^ H ♦ 2^ , W > ) - p » 



“ 2Di 2WJ_ - 2 j' tan JT - £_|_ + V (VJ 11 + 2V* £' ) 

f ;/-. ^ ir'v 2 



(4 v )r 



if • *?) | 



(71) 



D' 



D v v ' + % i ' ♦ D 5‘ " + ^>/° ' + % W' 



2D r 



‘ yt ♦ pi *1 

-V tf\ 



- 2 Di 



i - Pi 1 



j ~)v' + J( ' tan^ 

pL ♦ 



XI 



- W» » jp« 

^ f 



atmosphere parameter, P , and its derivative, /O » , may 

/ ' 



(72) 



The 

be evaluated as follows: 

For altitudes below 36,089 feet: 
/0 (1 - 6.89 x 10 -6 h )^256 



/”/ 



/ 



» - -It.256 /> Q (1 - 6.89 x 10“ 6 h )3.256 (6 C 89 x 10 -6 ^ 



r 



-asr 



But dh * tanY 
cbc 



Then: 



r 

/>» ■ - 29*38 tan^ 

y> 10 6 - 6.89h 



(73) 



/ 
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For altitudes above 36, 0 39 feet: 

* -f u , -sfcoaM 
y 5 - yO O (0„2971) e 



r ^ 



,(0.2971) 



-dh 

dx 

207^16 



And 






i > - tanY 

"2“7'S'6" 



(7b) 



'i’be formulation of the final equations for the unbounded 
case may now bo made as follows: 

Fy - F - F v ,» -Jy, ( tf'ta rtf - \ = 0 

V / 

Substuting for Fy, F/V, Fyi * , and Fyt from (55), (56), (57), 
and (53); and for Dy which occurs in Fy from (66) gives after rearranging: 



D 0 - Di 

T“ t* 



(3 -Si-/ 


y r 

- VcosJ 




g • 






tan$ 

(75) 



0 



From (35), tiie second Fuler equations: 

Fy + FtanV - F^\» - 1^^'taniS - V 1 j = 0 

Substituting for F^, Ftant, F yi* and Fjf from (59), ( 60 ) and (6l) 
and for the drag derivatives which are introduced from (69), (70) 
and (71) gives after rearranging: 

(D 0 - D^) tartf + '.J - 2 Dj I 2’ F< - . )(' tanX -/-* * 

oosi (-§ 7 *V) [ If 



7» + _V_ v >(Vy 11 + 27* S 1 ) J 
" s >♦*$'_ 



( 76 ) 






Solving (75) for V' and (76) for^" yields the two equations 



whose solution will define the optimum path* 

v ' “ - - g|~ - jiWi-Zgr)- 1 * 

wprv 2 cos^ ) \ 



I & 

i" 



(’-&) * 






vor 

g 



2Di 

- V* - y* tan)/ - p ' 

nr tp 



♦ 2g^' 


1 

tani 


J 

(' 




°T 




(D 0 - 


• Dj^) tani ♦ W 


- /It. 




cos)/ 


l V2 } 



- 2 V«tf« 

— F~ 



2W* 

nr 



(78) 



The final equations for the bounded case are formulated in 
like manner: 

From (1:3), the first Euler equation 



W“ - V* 



fy _ f y •“ Dy 

7 {£ - F-Wsin*/ 



- V" - 0 

nr 



Substituting for W n /W’ from (62) and (63) and for Dy and D' which 
then remain from (66) and (72) yields on rearranging: 

i. 



2 D, 



V' - 

V 



1 * PI 

V fjr 



2 Di 



F¥) 



- $ • ta.nl 



1 - V 1 



^ " * + Jil - A? 

+ V) w Yi 

V- ; ' 



- f* ♦ W'sin 






+ W^*cos)( * ( f' - v 'fy + W* \ (f - D — Wsin)() * 0 (79) 

If f W J 



From (1:6), the second Suler equation** 

Jjjl D^ 1 + (f - D - Wsin^) - f ^ - D^ - Wcos t ♦ D^i! 
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. ty, (f' - D' - W'siritf - w/'cosV ) -0 

* (f - ft - Wsin^) : 



Substituting for W"/W' from (62) and (63)? for Dy , i • ^ and D' 

from (70), (71) » and (72) yields on rearranging: 

f» + 3W' + V" - 2(f ' - W'sintf - VJ ^»cos^) - 2D' 

T" IT “V 1 ' (f - ft' - Wsin^j 



+ V2 _ (f - D-Wsin){) - f x + 2Di tan& - Wcos X 



♦ V(V)f" + 2V' /') - (?' - 2 ^'tan/ - 0 

' 4fV * £7 f 



(80) 



Solving (79) for " and (80) for V" yields two equations 
whose solution will define the desired path: 

'i n - (\ 1 1 If' - W'sW - W i'cosi - / f' - v ' f V + YJ» ) (f-D-Wsintf) 

ipssrj l v~ “7 w ) 



. -| v / 1 - 

- 2 Do [V±-l + p l)\ - 1 - VI) \ g J 

vr nr *7?/\ CT nr 77 + &7y£ t 

* • g / 



V" 



*- V'f* - 3W , V' + 2V'(f t -W'sin^ - W ^ * cos^ ) - 2V'D» 
~1 W (f - D - Wsin~6) 



(81) 



V' 



t / 

("fe + ^ ) I” f^ (f - D - Wsiriii)-- fy + 2B^tan/ - Wcos^ 

“TT5I [T * J 



- w« (V <" ♦ 2V' V«) + VV + 2V' ' tan/ 

(>•?) "T 



(82) 
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PART V: METHOD OF SOLUTION OF THE EQUATIONS FOR 

AN EXTREMAL PATH 

Considering first the unbounded or general equations for an 
aircraft flying a minimum fuel path* the following equations control 
the aircraft's motion: (21b), (65), (77) and (78)* These equations 
will be repeated here for ready reference: 

A. The general equation of motion along the flight path: 

F » D + W V' Vcosi * Wsintf 



g 

E« The drag equation: 



D - 1 PS V 2 ^ + 2(Wcos~5) 2 j 1 + ft) 2 

~Z~I ^eARyoS V 2 \ g / 



(A) 



(E) 



C. The first Euler equation: 



v ’ - £L 

worv 2 cos)i 




g r ♦ g°r + 2gV 

~V~j ~T 



D. The second Euler equation: 

■V* ■ 



')* 


(D 0 - &) tan# + W 






cos^ 

h» 


{ V2 ) L 



2W« 

W 



where 



- - {'tani - PI 

fj 



- 2V ' ^ • 

T 



(D) 



D o - 4-/> S v2 % 

■ 2(Wcos7f)2 ( 1 + A 1 yg )- 

'feAR/^S V 2 \ 8 ' 



This system is composed of four nonlinear differential equations 



Uo 



in o and V. Three of the equations are first order and one is second 
order. This is equivalent to a third order equation in one unknown 
and thus three initial or boundary conditions may be specified* 

Since no general closed form solution exists to this system, we 
are forced to integrate by means of a step-by-step process. This 
requires the use of initial conditions rather than boundary' conditions. 
From equations (A) and (B) it can be seen that the forces on the air- 
plane depend only on V, V’ , V , and 1 and therefore our initial 
values should be given in terms of three of these at a given altitude, h* 
Since V 1 can also be determined directly in terms of the other three 
from equation (C), it will be most convenient to specify V, ^ , and 
V 1 initially. 

For a given aircraft we know W 0 and Cp 0 , and for the unbounded 

case we assume unlimited thrust and knowledge of C? at any thrust* 

Now, if we specify V, 'i , and initially, V 1 is computed from (C) 
and D, D 0 and from (B). With V' and D known, F is determined 
from (A). Since C-j» is a function of F^ it was necessary to estimate 

initially to find V* . Now that F has been determined, an iterative 

procedure should be used to correct V', Gj> and F. V/* is found easily 

from: , 

VJ* ■ — C»p e (e) 

Vcost 

Now % M can be computed from equation (D). Thus we have defined all 
quantities for our starting point by initially specifying three 

variables V, ^ , and * *. 



By using the first two terms of a simple Taj lor series ex- 
pansion we may now integrate step-wise to point number 1 and so on in 



the following manner: 

Knowing V 0 ', ^o', W c ' and ^ 0 " we compute by Taylor series 




V 1 * V 0 ♦ V 0 'Ax 
i± ■ io + Ax 



ii' • i 0 id' Ax 



(F) 



Wi - W 0 + W 0 * Ax 



Also since 



h' « tan^ ' 

■ h Q * h' Ax 



(F) 



Now with Vy, ^l'» ^1> ^1 anc ^ equations (A) through (E) we can 



determine all quantities at point (l). Using Taylor series again we 



If greater accuracy is required than obtained by two term 
Taylor series, additional terms of the series may be used, the Milne method 
may be used, or the increment in x maj' be decreased. Decreasing 
the increment in x reduces the errors due to non-linearities of 
the function involved. The Milne method is described in reference 2 

In performing a step-bj-step solution to the unbounded case, 
several considerations should be kept in mind to insure that the path 
so determined will be a "flyable" ones 

1. ) The equations can possibly give a maximum fuel path instead 
of a minimum, 

2, ) Flight path angle, V , cannot be allowed to drive altitude to 
negative values or the a/c into stall. 



proceed to point (2) and so nn» 



3. ) Accelerations and decelerations must be within the physical 
capabilities of the aircraft* 

h.) Minor oscillations of the computed path in its initial 
stages are probably due to computational variables such as the 
size of x and to the initial conditions, and can be faired out for 
an actual flight path* 

Juggling of initial conditions and increment in x will help avoid 
these difficulties at the beginning of the patho 

To compute a minimum fuel path for the bounded case equations 
(C) and (D) must be replaced by (61) and (82). Equations (8l) and 
(82) will be restated here only in functional form: 

First Euler Equation 

f(V, i , V', V") - 0 (G) 

Second Euler Equation 

f(V, i , V', i*, V", ») - 0 (H) 

It must be remembered in the bounded case that thrust is 
fixed as a known function of V and ^ for a particular aircraft. 

With this in mind, a close imspection of all of the equations for the 

bounded case (i.e. (A), (B), (G), (H)) tells us that the second Euler 

equation (H) is not necessary to the solution of the system 0 This 
is as we would suspect from the physical conditions of the bounded 
problem, since, by setting thrust at a fixed rating, say maximum 
continuous, it becomes a known function of V andV. Thus one less 
equation is required for solving the sustem* 



h3 



This can be seen more clearly from a look at (G) and (ri), 
with the balance equations. We have two second order differential 
equations with second order terms on 9 in both (G) and (H) 
equations and second order terms on V in (G) only. y " could be 
eliminated between (G) and (H) and our system would be equivalent to 
a third order equation in one unknown. Thus three initial conditions 
may be specified. Looking at (G) alone with the balance equations 
we also have an equivalent third order equation in one variable. 

Thus three initial conditions may be specified in either case. 

If we initially specify V, y, and y' as before, for a particular 
aircraft, F, C^, and D are determined. Then since F and D are known, 
V* is determined from equation (A)© If V' is determined from 
equation (A), the V” term in equation (H) is redundant. The other 
variables in equation (G) are functions of known or specified 
variables so V, y, V', y‘ and y" are related uniquely in (G) and 
therefore (H) can be disregarded. Ey using only the first Euler 
equation all quantities are specified or determined for the initial 
point of the bounded problem. Py equations (F) representing the 
first two terms of the Taylor series, we can again integrate stepwise 

to point 1 and so on as before. 

, The bounded problem and also the more general unbounded problem 
each have specific applications which fit completely within the scope 
of one or the other of the two situations. It is felt, however, in 
the absense of any computed results that the general solution will 

CL 

be comoosite of the two. Due to the similarity of the systems of 
equations for the two cases, transition from one case to the other is 






an easy matter at any point oh the path. When thrust required ex- 
ceeds that available while using the unbounded equations, transition 
to the bounded system of equations should be made. This transition 
should be accomplished at the last point on the path at which the 

thrust requirement is equal to the thrust available* 

( 

Now, while proceeding step-by-step along the path computed by 
the bounded system of equations, it is necessary that each path point 
be checked in the unbounded eqriations. At that point on the path 
where the unbounded thrust requirement comes within the envelope of 
thrust, available from the power plant, transfer back to the unbounded 
system of equations should be made and computations continued until 
such time as the thrust requirement again forces transition to the 
bounded case. This procedure is continued until the required distance. 



x is coveredo 



CONCLUSIONS 



It is concluded that a system of equations for maximum range 
or minimum fuel flight paths may le written which are not restricted to 
level flight or to negligil le accelerations. The equations can be 
applied to flight paths which commence at take-off or at any other set 
of initial conditions from which an extremal path is desired. 

It is further conclude '- 1 that by specifying three initial 
conditions, methods of solution by stepwise integration exist for 
the system of equations. 

It is also concluded that, should the extremal path call for 
greater thrust than is availalle from the aircraft, a maximum range or 
minimum fuel path can be commuted based on maximum continuous thrust 
available. There are specific applications for loth types of paths 
individually as well as for a composite path. 

From Appendix A it is concluded that classical parameters for • 
maximum range are restrictive applications of the more general equations 



derived in this thesis<> 



APPEIimX A 



Reduction of Euler Equations to Classical Range Parameters 



Fraa the Euler equations representing a maximum range path at 
variable altitude it is a simple matter to derive equations or 
parameters for maximum range in level flight. Starting with equation (77) 
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Imposing the conditions for level flight ('/*'/' * 0) we have 
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Substituting for W 2 its level flight equivalent (W * -^sV^C^) 
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Assuming nominal values of Op, V, and C^, sa^: 



1 lb/lb hr * I/ 36 OO sec 
1:00 ft/sec 
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* in slightly different form 



(A.l) 



(A. 2) 



Then Oj V ■ *0035 and maj be considered negligible with respect to lo 

G 

And V'CjCtV 2 * *009 V' where V' is in ft/sec and is an extremely 

— F“ ft 

small quantity. This term na’- be considered to be zero. With these 
assumptions the equation reduces to: 



maximum range at any one given altitude. This parameter is perhaps not 
as familiar as the classical 
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Or: 
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Equation (A #3') represents a parameter defining the speed schedule for 
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parameter but it can be shown that the two are equivalent* 
Differentiating " K 
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from which 
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APPENDIX B 



Consideration of Mach Number Effects 

In the foregoing developement of maximum range and minimum 
fuel equations, only the sub-sonic speed region was considered. In 
this section, Mach effects on the equations will be determined and 
applied, thereby removing any speed limitations, 

Mach number appears explicitly in the drag terms of both the 
unbounded and bounded cases, and may possibly appear in the thrust, f, 
of the bounded case. Therefore, by determining the Mach corrections 
to the D and f terms and their derivatives, the complete system of 
equations is corrected and is applicable to the entire speed range. 

Considering drag first, the Mach effects thereon may be accounted 
for as follows: 
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Now Cp o and e are functions of Mach number. Defining 
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Drag mA'v written as: 
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Taking derivatives on drag: 
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Since neither drag nor Kach number are functions of V', as before: 

. “ 0 also since * c>M “0 

^ -sr 

the derivatives on drag with respect to ^ and are unchanged 
from their subsonic form provided Cp and e are inserted as proper 




functions of Mach number* 

To evaluate Py' , from (70) and the above discussion: 
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From (71) and (E2) 
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To evaluate D* 
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Then from (gl) and (72) 
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As was previously noted, the bounded thrust, f, may possibly 
be a function of Mach number. If this is the case, derivatives 
may be taken in the normal manner as follows: 
f - f(V, h, K) 
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f “ 6t c)M V* ♦of V» ♦ 3f dh x ' 
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The partial derivatives with respect to the various variables being 
obtained from known power plant characteristics© 
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